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Transport of dipolar Bose-Einstein condensates in a one- dimensional optical lattice 
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We show that magnetic dipolar interactions can stabilize superfluidity in atomic gases but the 
dipole alignment direction required to achieve this varies, depending on whether the flow is oscillatory 
or continuous. If a condensate is made to oscillate through a lattice, damping of the oscillations can 
be reduced by aligning the dipoles perpendicular to the direction of motion. However, if a lattice is 
driven continuously through the condensate, superfluid behavior is best preserved when the dipoles 
are aligned parallel to the direction of motion. We explain these results in terms of the formation 
of topological excitations and tunnel barrier heights between lattice sites. 
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I. INTRODUCTION 

Since the development of optical lattices for cold atoms 
[l| there has been a wealth of research in this area, not 
least on transport properties. Recent developments in 
the cooling of atoms with strong magnetic dipolar inter- 
actions |2|-l5|] have broadened perspectives further, with 
the possibility to explore how magnetic effects influence 
superflow; analogous problems in solid state physics have 
been a matter of debate for many years [6] . 

To date there have been a number of experimental 
studies on continuous and oscillatory transport of cold 
atoms in periodic potentials with non-magnetic atoms 
[7] , including work on Bloch oscillations [8, 9] , critical 
velocities and dynamical instabilities [10|, [ll| • The first 
experiments with dipolar gases in optical lattices are 
now starting to reveal information about static properties 
[l2[ . Theoretical work on gases in optical lattices has pri- 
marily focused on static phase diagrams pj| and dipolar 
interactions have also been considered in such cases [1^ . 
Dynamical effects in trapped dipolar gases without a lat- 
tice potential have been studied [15], such as excitation 
spectra [l^, ^^ , and anisotropic superfluidity [18|, [l9| . 

However, there has been little work on the transport of 
dipolar gases in optical lattices [20]. In particular there 
has been no consideration of how the directional align- 
ment of magnetic dipoles in such systems affects the su- 
perfluid properties. There has been theoretical work on 
lattice transport employing band structure calculations 
in one-dimension [211, |22| but these are difficult to ex- 
tend to many experimentally relevant situations in higher 
dimensions, especially when non-local interactions must 
be taken into account [231]. In addition such approaches 
leave out certain details of the dynamical evolution of 
the cloud [24]. Saito et al. used a Hubbard model with 
a Gutzwiller Ansatz to show that dipoles aligned per- 
pendicular to a two-dimensional lattice plane had only 
a small effect on dynamical instability thresholds, but 
other configurations were not considered. 

Here we show, using dynamical simulations of the 
non-local Gross-Pitaevskii equation, that the directional 



alignment of magnetic dipoles in a Bose-Einstein conden- 
sate plays a crucial role in determining its transport prop- 
erties through a shallow one-dimensional optical lattice. 
Interestingly, the direction required to best stabilize the 
superfluidity varies depending on whether we consider 
oscillatory flow or continuous flow. We show this occurs 
because the physical processes that break down the su- 
perflow differ in the two cases. In the case of oscillatory 
flow, superfluidity is broken down by the formation of 
topological excitations. In the case of continuous flow no 
such excitations are observed, but dipolar interactions 
act to increase or decrease the tunneling barrier between 
lattice sites, depending on the dipole orientation. 

The paper is organized as follows: in section |II] we ex- 
plain the numerical methods to solve the non-local Gross- 
Pitaevskii equation for a dipolar gas in an optical lattice. 
In section [nil we present our results for continuous and 
oscillatory flow. Finally, we conclude in section HVl 



II. NUMERICAL APPROACH 

In both the oscillatory and continuous cases we model 
the dynamics of the dipolar Bose-Einstein condensate 
(BEG) using the full three-dimensional Gross-Pitaevskii 
equation (GPE) [25|, [2y], including the usual integral 
term for the dipolar potential |27j 



^^^V^(r,t) 
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where '^(r, t) is the wave function at position r and time 
t, normalized to the total particle number N = 10^, ra 
is the mass of a single particle and g = Anh'^a/m is the 
contact interaction strength with a the 5- wave scattering 
length; the other symbols have their usual meaning. The 



external potential 

Kxt(r) = l^osin2(^) 



+ -m (uolx^ + 



2 2 



■^Iz') 



(2) 



includes the optical lattice with spatial period d and 
depth Vq^ and the harmonic trapping potential with the 
trapping frequencies uox^ ^y and uoz- The potential Vdd{^) 
describes the interaction between two polarized magnetic 
dipoles and is given by 



Kid(r) = 



/io|/L6pl-3cos^(i9) 

47r IrP 



(3) 



where 6 is the angle between an external magnetic field 
which polarizes the dipoles and the vector r between the 
position of the two dipoles (we assume perfect alignment 
with the field), ii is the dipole moment and /io is the 
vacuum permeability. We choose ^^Cr atoms for our 
simulations, as in recent experiments [12] which have 
ImI = 6/iB where /iB is the Bohr magneton and we set 
a = 2.0 X 10~^ m [28]. For the optical lattice we assume 
d = 1.59 /im and Vq = 1.53 Er {Er = Ti^t:'^ /2md^ is the 
recoil energy) . The frequencies of the harmonic trap are 
chosen to be {ujx^^y^^z) = 27r x (48 Hz, 16 Hz, 16 Hz). 
With these parameters we obtain a system that is deep 
in the superfiuid regime. 

To solve Eq. [TJ we use a standard Fourier split-step 
method where we treat the dipolar term with the convo- 
lution theorem [23, [29| . We prepare the numerical BEC 
ground state by starting with the Gaussian analytical 
solution for the harmonic trap [30] and ramping up the 
interactions and the lattice potential adiabatically. 

For each model system we look at three configurations 
to determine the infiuence of the dipole-dipole interaction 
(DDI). For the first configuration we neglect the DDI and 
look at the system with pure contact interactions which 
will be referred to as the Contact Configuration (this sys- 
tem is fictitious and is used purely for academic compar- 
ison). Secondly, we look at the system with the dipoles 
polarized along the x-direction (parallel to the direction 
of motion), which we call the Parallel Configuration. Due 
to the geometry of our model system the DDI raises the 
total potential energy in this case. Thirdly, we look at the 
system with the dipoles polarized along the ^/-direction 
(perpendicular to the direction of motion) , henceforth re- 
ferred to as the Perpendicular Configuration. In this case 
the DDI reduces the total potential energy. In the follow- 
ing section |TTl] we always depict the results for the Con- 
tact Configuration with red solid curves, and the results 
for the Parallel Configuration (Perpendicular Configura- 
tion) with green dashed curves (blue dotted curves). 



III. RESULTS 

Here we present our results for oscillatory fiow and con- 
tinuous fiow. We are primarily interested in the break- 
down (or preservation) of superfiuidity depending on the 




FIG. 1. (Color online) Sketch of the model system for oscil- 
latory flow. The trap is displaced a distance Ax in the x- 
direction to trigger a center-of-mass oscillation. The colored 
oval represents the BEC, the gray line shows the trap before 
the displacement, the black line the trap after displacement 
and the oscillating dark red line the optical lattice. Black 
arrows indicate the axes with r the radial direction. 



direction of the magnetic dipoles. Unfortunately, the su- 
perfiuid fraction in our systems cannot be unambiguously 
defined in a useful, quantitative way. This is primarily 
due to the non-equilibrium and inhomogeneous nature 
of the system J3lM36|. Naive use of the Landau crite- 
rion is further complicated by the spatially extended na- 
ture of the lattice potential that perturbs the system [30| . 
We therefore make use of more qualitative measures. In 
the oscillatory case, we consider the damping of in-trap 
dipole oscillations as has been done previously [37]. In 
the continuous case, we consider the total energy of the 
BEC; a superfiuid gas will acquire no energy from the 
moving lattice, therefore an increase in the BEC energy 
indicates a breakdown of superfiuidity. 



A. Oscillatory flo^v 

To create oscillatory flow we displace the trap a dis- 
tance Ax along the x-direction at t = 0.0 and trigger a 
center-of-mass oscillation as sketched in Fig. [TJ Due to 
the harmonic trap the BEC oscillates back and forth with 
a velocity that varies sinusoidally (we do not fulfill the 
conditions for Bloch oscillations, a displacement of ap- 
proximately 7 /im would be required for the condensate 
to reach the edge of the first Brillouin zone). In Fig. [2] 
we show the expectation value of the BEC's position in 
the x-direction as a function of time for a trap displace- 
ment of Ax = 4.5 /im (this gives a velocity amplitude of 
cJa^Ax ~ 1.4mms~^). The Contact Configuration (red 
solid curve) shows slight damping of the center-of-mass 
motion during the simulation. However, the Perpendicu- 
lar Configuration (blue dotted curve) shows less damping 
compared with the red solid curve. By contrast, the Par- 
allel Configuration (green dashed curve) shows stronger 
damping. To further investigate the system we consider 
snapshots of the density profile following several oscilla- 
tions in the trap. Figure [3] shows the disruption of the 
BEC density profile as a result of oscillatory fiow. Panel 
(0) shows the profile at the start of the simulation before 
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FIG. 2. (Color online) Expectation value of the BEC position 
along the x-direction for the three configurations: without 
dipolar interactions (red solid curve), with dipoles polarized 
parallel to the direction of motion (green dashed curve) and 
with dipoles polarized perpendicular to the direction of mo- 
tion (blue dotted curve). 




FIG. 3. (Color online) Cross-section through the z — ^ plane 
of the atom density distribution after t — 0.08s oscillatory flow 
for the configuration without dipolar interactions (a), with 
the dipoles polarized parallel to the direction of motion (b), 
and with the dipoles aligned perpendicular to the direction of 
motion (c). For comparison we also show the profile at t = 0.0 
without dipolar interactions (0). 



flow has begun (t = 0.0). We use the Contact Configura- 
tion here but differences due to the influence of the DDI 
are barely discernible at this point. Panels (a-c) show 
the density profile after 0.08 s oscillatory flow (around 4 
oscillation periods). Panel (a) shows the Contact Config- 
uration; we see that a limited amount of disruption has 
occurred, compared with the situation at t = 0.0, the 
cloud acquiring a small "tail" to the left. By contrast, 
panel (b) shows the Parallel Configuration. We see sig- 
nificantly more disruption in this case, including a large 




FIG. 4. (Color online) Sketch of the model system for contin- 
uous flow. The colored oval represents the BEC. The optical 
lattice is dragged along the x-axis through the cloud at con- 
stant velocity. The black line shows the harmonic trap, the 
oscillating gray curve shows the optical lattice at t = 0.0 and 
the oscillating dark red curve shows the shifted lattice. Black 
arrows indicate the axes with r the radial direction. 



density minimum in the center, and radial topological 
excitations. These topological excitations typically con- 
sist of solitons and vortices and form as a result of the 
direction reversal of the BEC's motion [38] which pushes 
atoms into the radial direction r; this nonlinear mixing 
explains the damping behavior in Fig. [2l Panel (c) shows 
the profile for the Perpendicular Configuration. In this 
case, the "tail" is smaller and the cloud has largely re- 
tained its form, compared with the starting profile. 

The reason we see less damping for the Perpendicular 
Configuration is because the DDI creates potential wells 
that restrict motion in the radial direction and therefore 
the cloud is more stable against motion away from the 
primary direction of flow. By contrast the Parallel Con- 
figuration amplifies the creation of topological excitations 
because the DDI raises the potential energy, and thereby 
the level of instability. 



B. Continuous floAv 

In our second system we drag the optical lattice 
through the BEC at constant velocity v, as sketched in 
Fig. m Figure [5] shows the total energy per atom as 
a function of time for different v. We see that below 
V ~ 1.00 mm s~^, little energy is transferred to the BEC 
[Panel (a)]. However, around this velocity, the three con- 
figurations begin to gain energy [Panel (b)]. Initially, this 
energy gain is similar in all three cases but for higher 
velocities, the curves start to diverge from one another 
[Panel (c) and Panel (d)]. The Perpendicular Configu- 
ration (blue dotted curve) shows a greater increase in 
energy than the others whereas the Parallel Configura- 
tion (green dashed curve) has the smallest energy in- 
crease. The Contact Configuration (red solid curve) lies 
between the other two configurations. We see that the in- 
crease in energy only becomes significant after ~ 10 ms. 
For this reason, these instability effects are unlikely to 
have played a significant role in the oscillatory case, even 
though the velocities were above the required thresholds; 
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FIG. 5. (Color online) Total energy per atom as a function 
of time for lattice velocities 0.75 mm s~^ (a), 0.98 mm s~^ (b), 
1.06 mm s~^ (c) and 1.13 mm s~^ (d). The red solid curves 
show the configuration without dipolar interactions, the green 
dashed curves the configuration with the dipoles aligned par- 
allel to the direction of motion and the blue dotted curves the 
configuration with the dipoles aligned perpendicular to the 
direction of motion. 



the condensate spent too little time above the thresholds 
for the atoms to respond. 

As in the case of oscillatory flow, we explore this be- 
havior further by examining the cloud density profiles. 
Figure [6] shows the density profiles of the gas after 0.04 s 
flow time with v = 1.13 mms~^. Panel (0) shows the 
starting profile at t = 0.0 as before. Panel (a) shows 
the profile for the Contact Configuration. We see that 
significant numbers of atoms have been displaced from 
their original positions in the center. We are at the edge 
of the zone for dynamical instability where perturbations 
in both the -\-x and —x directions can grow exponentially 
in time. The density of these "wings" is noticeably less 
for the Parallel Configuration [Panel (b)]. By contrast, 
the Perpendicular Configuration shows increased density 
in the wings [Panel (c)]. Unlike the oscillatory flow, we 
see no large density minima and little evidence of radial 
topological excitations. 

To summarize these results, we show the total energy 
per atom as a function of the lattice velocity in Fig. [3 
The values were taken after t = 0.02 s which is well above 
the correlation time of the BEC. The energy increase is 
noticeably more severe for the Perpendicular Configura- 
tion than the Parallel Configuration. For comparison, in- 
stability thresholds have been calculated by determining 
the lowest energy band and checking its stability against 
perturbations following reference [22] with pure contact 
interactions [39]. The thresholds for energetic and dy- 
namical instability from the band structure calculations 
correspond well to all three of our configurations as ex- 
pected [33 • 

That the Perpendicular Configuration can stabilize su- 
perflow for oscillatory motion but destabilize it for con- 
tinuous flow (vice-versa for the Parallel Configuration) 
may initially seem paradoxical. However, we have ob- 
served that the superflow is not broken down by topolog- 




FIG. 6. (Color online) Cross-section through the z = plane 
of the atom density distribution after t = 0.04 s continuous 
flow for the configuration without dipolar interactions (a), 
with the dipoles polarized parallel to the direction of motion 
(b) and with the dipoles aligned perpendicular to the direction 
of motion (c). The lattice velocity is v = 1.13 mms~^. For 
comparison we also show the profile at t = 0.0 without dipolar 
interactions (0). 



ical excitations in the continuous case, as it is in the oscil- 
latory case. We therefore turn our attention to tunneling 
between lattice sites by considering the effective lattice 
potential 14ff(r) which is the sum of the optical lattice, 
the contact interaction and the dipolar interaction. If we 
compare the depth of Kff (r) at x = for the Perpendic- 
ular Configuration with the Contact Configuration, we 
find that the DDI increases the effective lattice depth by 
about 3.0%, making it harder for atoms to tunnel. Con- 
versely, the Parallel Configuration lowers the depth of 
the effective potential by about 3.6%, effectively increas- 
ing the hopping coefficient. Consequently, the reduction 
(increase) of the tunnel barrier facilitates (hinders) the 
motion of the BEC through the lattice, thereby explain- 
ing the results in Figs. [5]- [3 



IV. CONCLUSIONS 

In conclusion, we have studied superfluid transport 
of magnetic dipolar Bose-Einstein condensates in one- 
dimensional optical lattices. We have shown that the 
dipolar interactions can be used to stabilize the superflow 
in such systems but the atomic dipoles must be aligned in 
different directions to achieve this, depending on whether 
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FIG. 7. (Color online) Total energy per atom at t = 0.02 s 
as a function of v for the three configurations: without dipo- 
lar interactions (red dots), with the dipoles polarized par- 
allel to the direction of motion (green diamonds) and with 
the dipoles polarized perpendicular to the direction of mo- 
tion (blue squares). Vertical gray lines show the thresholds 
for the energetic instability (dashed curve) and the dynamical 
instability (dot-dashed cure). 



the flow of the gas is oscillatory ( "alternating current" ) or 
continuous ( "direct current" ) . In spite of the apparently 
paradoxical nature of this result, we showed it could be 
simply explained by the dipoles' effect on topological ex- 
citations and optical lattice barrier height. The direction 
reversal in the alternating case is the key difference that 
gives rise to different physics. 

It is likely that the effects we studied here will be even 
more pronounced for cold atoms with stronger magnetic 
dipole moments such as Erbium and Dysprosium, assum- 
ing that possible instabilities can be controlled. There are 
further questions relating to two- and three-dimensional 
lattices; interesting phenomena may occur in deeper lat- 
tices. In addition, there is the possibility of considering 
finite-temperature gases and studying the effects of mag- 
netic dipoles on the superfluid critical temperature. 

We thank T. Pfau and A. Griesmaier for helpful dis- 
cussions and cooperations within the SFB TR21. We 
also thank R. Scott and A. M. Martin for discussions, 
and BW-Grid computing resources. 



[1] 
[2] 
[3] 

[4] 
[5] 

[6] 
[7] 

[9: 
[lo: 

[11 

[12 

[13 

[14 
[15 



O. 



Morsch 



and 



M. 



Oberthaler, 



Rev. Mod. Phys. 78, 179 (2006) , 

A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and 

T. Pfau, Phys. Rev. Lett. 94, 160401 (2005)| 

Q. Beaufils, R. Chicireanu, T. Zanon, B. Laburthe-Tolra, 

E. Marechal, L. Vernac, J.-C. Keller, and O. Gorceix, 



Phys. Rev. A 77, 061601 (2008) 



M. L u, N. Q. Burdick, S. H. You n, and B. L. Lev, 

Phys. Rev. Lett. 107, 1904 01 (2011)[ 

K. Aikawa, A. Frisch, M. Mark, S. Baier, 

A. Rietzler, R. Grimm, and F. Ferlaino, 
|Phys. Rev. Lett. 108, 210401 (2012) . 

W. Zhu, H. Kang, Y. C. L ee, and J. C. Chen, 
"Phys. Rev. B 50, 10302 (1994): 

J. H. Denschlag, J. E. Simsarian, H. Haffner, C. McKen- 
zie, A. Browaeys, D. Cho, K. Helmerson, S. L. Rolston, 
and W. D. Phillips, J. Phys. B^35, 3095 (2002), 

B. P. Anderson and M. A. Kasevich, 

Science 282, 1686 (1998)] 

O. Morsch, J. H. Miiller , M. Cristiani, D. Ciampi ni, and 
E. Arimondo, Phys. Rev. Lett. 87, 140402 (2001) , 
S. Burger, F. S. Cataliotti, C. Fort, F. Minardi, 
M. Inguscio, M. L. Ch iofalo, and M. P. Tosi, 
|Phys. Rev. Lett. 86, 444r(2001 ) . 

L. Fallani, L. De Sarlo, J. E. Lye, M. Mod- 
ugno, R. Saers, C . Fort, and M. Inguscio, 

|Phys. Rev. Lett. 93, 14 0406 (2 004) , 
S. Miiller, J. Billy, E. A. L. Henn, H. Kadau, A. Gries- 
maier, M. Jon a-Lasinio, L. S antos, and T. Pfau, 
[Ph ys. RevTA 84, 053601 (201 1)[ 

D. v an Cost en, P. van der St rat e n, and H. T. C. St oof, 
Phys. Rev. A 63, 053601 (2001)f 

K. Goral, L. Santos, and M. Lewenstein, 

Phys. Re v. Lett. 88, 170406 (2 002) , 

R. M. W. van Bijnen, N: G. Parker, S. J. J. M. F. 
Kokkelmans, A. M. Martin, and D. H. J. O'Dell, 
[Phys. Rev. A 82, 033612 (2010) [ 



[16 

[17: 
[is: 

[19 

[2o: 

[21 

[22: 

[23: 

[24 

[25: 
[26: 



[27 

[28: 
[29: 



[30 



[31 

[32: 



L. Santos, G. V. Shlyapnikov, and M. Lewenstein, 

Phys. Rev. Lett. 90, 250403 (2003) , 

S. Ronen, D. C. E. Bortolotti, and J. L. Bohn, 

fPhys. Rev. A 74, 01 3623 (2006) , 

R. M. Wilson, S. Ronen, a nd J. L. Bohn, 

[Phys. Rev. Lett. 104, 094501 (2010)] 

C. Ticknor, R . M. Wilson, and J. L. Bohn, 

| Phys. Rev. Letl 106, 065301 (2011), 

M. Fattori, G. Roati, B. Deissler, C. D'Errico, M. Za- 

ccanti, M. Jona-Lasinio, L. Santos, M. Inguscio, and 

G. Modugno, [Phys. Rev. Lett. 1017190405 ];2068j[ 



B. 



Wu, R. 



B. 



Diener, 



and Q. 



Niu, 



[Phys. Rev. A 65, 025601 (2002) , 
M. Machholm, C. J. Pethick, 
Phys. Rev. A 67, 053613 (2003) , 
Y. Lin, R.-K. Lee, Y.-M. Kao, 



and H. Smith, 



and T.-F. Jiang, 



[Phys. Rev. A 78, 023629 (2008), 

V. Brazhnyi and V. Konotop, 

[Mod. Phys. Lett. B 18, 627 (2004 ). 
T. E. Judd, R. G. Scott, and T. M. Fromhold, 
Phys. Rev. A 78, 053623 (2008): 

T. E. Judd, R. G. Scott, G. Sinuco, T. W. A. Mont- 
gomery, A. M. Martin, P. Kriiger, and T. M. Fromhold, 

[New J. Phys. 12, 063033 (2010)> 

K. Goral and L. Santos,|P hys. Rev. A 66, 023613 (2002)[ 
The scattering length can be tuned through a wide range 
using Feshbach resonances. This is a typical value [4Q| |. 
T. Lahaye, J. Metz, B. Frohlich, T. Koch, M. Meister, 
A. Griesmaier, T. Pfau, H . Saito, Y. Kawaguchi, and 
M. Ueda, |Phys. Rev. Lett. 101, 080401 (2008) [ 

C. J. Pethick and H. Smith, 

\Bose- Einstein condensation in dilute gases, 2nd ed. 
(Cambridge Univ. Press, Cambridge, 2008). 
O. Penrose and L. Onsager, Phys. Rev. 104, 57 6 (1956)1 
K. Goral, M. Gajda, and K. Rzazewski, 

[Phys. Rev. A 66, 051602 (2002) 



[33] R. Dodd, C. Clark, M. Edwards, and K. Burnett 



|Opt. Express 10, 284 (1997)' 
[34] P. B. Blakie and M. 

|Phys. Rev. A 72, 063608 (2005)] 
[35] K. Huang and H.-F. 

|Phys. Rev. Lett/eQ, 644 (19927| 
[36] C. Krumnow and A. 

Phys. Rev. A 84, 02160 8 (2011), 
[37] T. Saito, I. Danshita , T. Oz aki 

|Phys. ReVTA 86, 023623 (2012) 1 
[38] R. G. Scott, A. M. Martin, S. Bujkiewicz, T. M. 



Davis, 
Meng, 
Pelster, 
and T. Nikuni, 



Fromhold, N. Malossi, O. Morsch, M. Crist iani, and 
E. Arimondo, Phys. Rev. A 69, 033605 (2004) . 

[39] We have not performed these calculations including the 
dipolar interactions for the numerical reasons already dis- 
cussed. 

[40] T. Lahaye, T. Koch, B. Frohhch, M. Fattori, 
J. Metz, A. Griesmaier, S. Giovanazzi, and T. Pfau, 
[Nature 448, 672 (2007)] 



